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Abstract 



We made a classical trajectory Monte Carlo (CTMC) calculation of state 
selective cross sections for processes between some light ions and excited 
helium. The results, useful for analysis of spectroscopic data of fusion devices, 
are in good agreement with theoretical predictions of scaling laws. 



1 



To be published in Phys. Rev. A 



2 



Impurities in high temperature plasma devices are mostly found as fully 
stripped ions. No spectral line is emitted by these ions, but they may be 
detected through characteristic lines emitted by one-electron ions found in 
excited states after a process of charge exchange (CX) with a neutral atom: 

N + X''^N+ + X''-\n') ^N+ + X''-\n) + hv (1) 

where N is the neutral atom, X is the ion, n and n' are quantum numbers, 
and hv is the emitted photon. 

In fusion devices beams of hydrogen, helium and lithium are usually in- 
jected into the plasma to heat it or to monitor such plasma parameters as ion 
density, temperature, velocity of rotation [1]. Extraction of useful informa- 
tions from emission measurements relies on knowledge of the effective cross 
section for process (1). 

Process (1) may take place with in the ground or in an excited state: 
charge exchange involving neutrals in the ground state predominantly pop- 
ulate hydrogen-like states with principal quantum number n < 5 for low-Z 
ions, up to oxygen [1]. The decay of these states occurs mostly through 
An = 1 transitions, corresponding to emission in the X-ray or ultraviolet 
region [1]. The emission in the visible region is originated from states with 
higher n (e.g. transitions 10 ^ 9 at 6070 A of 0^+ or 8 ^ 7 at 5292 A of 
C^+); in this spectral range the detection of radiation is simpler and multi- 
chord systems, to spatially resolve the emission, are usually working. Cross 
sections for electron capture into these states from ground state neutrals are 
much smaller, but, on the other hand, the contribution from excited neutrals 
is greatly enhanced [1] so, even if the atoms excited by Coulomb collisions 
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with plasma particles constitute a small fraction of the entire beam, they 
may be the major responsible for the emission. This implies that to deduce 
informations from CX emission it is necessary to have data on the processes 
involving excited neutrals. 

In this work wc propose to calculate partial cross sections for process (1) 
(because of their much smaller cross section double electron processes are not 
here considered) between fully stripped ions and helium in the first excited 
state. We focus on conditions frequently encountered in plasma devices : 
the collision energy is varied in the range below 110 KeV/amu and the ionic 
charge varies from 3 to 8. 

A fully quantum mechanical treatment of the scattering process for such 
a system would be difficult because of the large number of states involved. 
The classical trajectory Monte Carlo (CTMC) method , where the quantum 
system is replaced by a classical one with all the particles interacting through 
Coulomb forces and their motion as governed by Newton's law, has proven 
to be a feasible and reliable tool in this kind of calculations. 

The system helium+ion actually represents a four-body problem, but the 
excited electron of helium plays a major role, so we chose to substitute the 
other electron and the nucleus altogether with a point- like particle interacting 
with the other charges through a central model potential [2] , which in atomic 
units (a.u.) can be written: 

V„„,W = -i-(i±^^H!:^exp(-2.499r). (2) 

r r 

The first eigenvalues of the Schrodinger equation with this potential are 
—0.904 a.u. and —0.156 a.u.. It may be proved, following similar lines 
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as in [3], that Vmod provides a good approximation of electron moment dis- 
tribution, which is a critical requisite for every CTMC computation [3,4] . 

Initialization of the system was performed according to the procedure 
developed by Reinhold and Falcon [3]. The initial distance between nuclei 
was set to lOZ, but checks were made to larger distances, up to 20Z, to 
ensure independence of the results from this parameter. For each energy a 
number of runs ranging from 4 • 10^ to 27 • 10^ was made to have a reasonable 
statistical error (< 10 per cent for the processes characterized by the greater 
frequencies) . We have considered ions which are of interest in fusion devices: 
carbon and oxygen usually constitute the most abundant impurities, while 
beryllium and boron are often used to condition vacuum chamber inner walls. 
Only distribution on the principal quantum number n was considered: at 
the densities and quantum numbers considered there is a complete Z-mixing 
before spontaneous decay take place [5] . 

In Fig. 1 are plotted electron capture cross sections vs n for different 
energies. The position of the maximum in the n~distribution of the cross 
section closely follows the rule: 

^3/4 

where / = —0.156 a.u. is the ionization potential of the helium. Relation (3) 
is already known to hold for collisions between ions and H or Li beams. 

If compared with analogue curves plotted in [6] relatively to collisions 
between ions and He in the ground state, we notice that our distributions 
are much broader. The spreading in n is due the increasing of the ratio Vp/ve, 
where Vp is the impact velocity of ion, and Ve is the initial orbital velocity of 
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the electron [5]. 

Scaling formulas with respect to charge of the ion, principal quantum 
number of the neutral, and relative energy of incidence have been derived for 
electron capture processes involving hydrogen. We verified all our data lie 

on a single curve after the rescaling 

(see Fig. 2). Scaling (4) has been suggested, for instance, in [7]. Our scaling 
well resembles those obtained in similar calculations with H or He in the 
ground state (for a comparison see the curves plotted in [6]). In our case n 
in eq. (4) is the effective (noninteger) quantum number. 

In alternative to process (1) single ionization is also a possible result: 

AT + X« ^ + X« + e (5) 

For ionization a similar relation to (4) has been suggested in [8] to hold, as 
a result of a close-coupling calculation on ground state and excited lithium: 

<'UE)--^, E^^, (6) 

Even though our plots cannot discriminate little differences in the exponent 
of Z, we may state that our results are very similar to those of [8]; our 
reduced cross sections are plotted in Fig. 3. 

In order to make sure that Vmod i^ay represent a good enough potential 
we repeated some of the calculations by using a more sophisticated potential: 
in Fig. 4 we plotted partial cross sections for the scattering with Beryllium at 
10 KeV/amu as calculated with Vmod and with the model potential reported 
in references [9,10], the agreement being quite good. 
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Summarizing, we have extended a well known method of calculation to 
processes not previously considered, which are of importance for the diagnos- 
tics of high temperature plasmas. Our results are in good accordance with 
extrapolations from previous calculations and empirical laws, so we may be 
confident on their accuracy. 

This work has been partially supported by the Ministero dell'Universita 
e della Ricerca Scientifica e Tecnologica (MURST). We wish to thank Dr. 
M.E. Puiatti and Dr. L. Carraro for their careful reading of the manuscript. 
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Figure Captions 



Figure 1: n-distribution for one electron capture between helium and beryl- 
lium, boron, carbon and oxygen. Cross sections are in 10^^^ cm^. Dotted 
curve correspond to an impact energy of £^ = 30 KeV/amu, barred curve to 
E — 20 KeV/amu and solid curve to = 10 KeV/amu. 
Figure 2: Reduced cross section a/ [n^Z] for one electron capture vs reduced 
energy Er? jZ^I'^ . Cross sections are in 10~^^ cm^. 

Figure 3: Reduced cross section o"/(n^Z^'^) for single ionization reduced 
energy Er? jZ^I"^ . Cross sections are in 10~^^ cm^. 

Figure 4: Partial cross sections a for single electron capture from Beryl- 
lium at 10 KeV/amu obtained through use of two different model potentials. 
Dashed line: Vmod of equation (2); solid line: model potential from [9,10]. 
Cross sections are in 10~^^ cm^. 
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